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Abstract. We introduce a new Macaulay2 package, Nauty, which gives access to powerful 
methods on graphs provided by the software nauty by Brendan McKay. The primary 
motivation for accessing nauty is to determine if two graphs are isomorphic. We also 
implement methods to generate families of graphs restricted in various ways using tools 
provided with the software nauty. 



Introduction. Let G and H be two finite, simple, undirected graphs on the common vertex 
set V with edge sets E(G) and E(H), respectively. We say that G and H are isomorphic 
if there is a bijection ip from V to itself which preserves edges, that is, {it, v} G E(G) if 
and only if {ip(u), <p(v)} G E(H). Determining whether two given graphs are isomorphic is 
known as the Graph Isomorphism problem. 

When Garey and Johnson wrote their classic book [GJ] on the complexity of algorithms, 
they specified twelve problems of ambiguous complexity, one of which was the Graph Iso- 
morphism problem. Unfortunately, it is still unknown if the Graph Isomorphism problem is 
P or NP-complete. Moreover, the problem is of such notoriety that some have even begun 
referring to a new complexity class, GI, of problems which reduce in polynomial time to 
the Graph Isomorphism problem [J]. Despite this, there exists computer software which is 
capable of determining whether two graphs are isomorphic in reasonable time. One such 
piece of software is nauty [N] by McKay. 

The nauty software is written in highly portable C and is designed to, above all else, 
compute whether two graphs are isomorphic. It also includes an extensive family of tools, 
collectively called gtools, to generate systematic modifications of graphs, to generate specific 
families of graphs, to generate random graphs, to filter a set of graphs for given properties, 
and to canonically relabel graphs. Most, if not all, of these features would be beneficial to 
any computer software that handles graphs. 

The package Edgeldeals [FHT] by Francisco, Hoefel, and Van Tuyl implements structures 
and methods for manipulating graphs (and hypergraphs) within Macaulay2 [M2], a software 
system by Grayson and Stillman designed to aid in research of commutative algebra and 
algebraic geometry. We introduce a new package, Nauty, for Macaulay2, which provides an 
interface with nauty. 1 Most of the aforementioned tools in gtools are accessible through 
Nauty. In particular, the methods of perhaps the greatest interest are arelsomorphic, 
f ilterGraphs, generateGraphs, and generateRandomGraphs. 

The remainder of this note is broken in to two sections: the first describes briefly the 
theoretical underpinnings of nauty and the second gives an example session of using Nauty 
along with a few useful caveats. 
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Canonical labellings. In [M], McKay describes the improved algorithms which he de- 
veloped to canonically label a graph; these algorithms are the heart of nauty and are sum- 
marised in [K]. We recall briefly the theoretical ideas which make such algorithms useful. 

Let G be a finite, simple, undirected graph on the vertex set V . We call an ordered 
partition ir = (Vi, . . . , V m ) of V a colouring of G and call (G, n) a coloured graph. Given 
a permutation a of V, we define ciyi) to be the set {o~{y) \ v G Vi}, a(n) to be the 
colouring (a(Vi), . . . ,a(V m )), and o~(G) to be the graph on V with edge set E(a(G)) given 
by {a(u),a(v)) G E(a(G)} for all {u,v} G E{G). If o~(ir) = ir, then a is called colour- 
preserving. 

If V = [n] = {1, . . . , n} and 7r = (Vi, . . . , V m ) is a colouring of V, then c(7r) is the colouring 
({1, . . . , |Vi|}, {|Vi| + 1, • - • , |Vi| + IV2I}, - • - , {n— \V m \ + 1, . . . , n}). A canonical labelling map 
is a function I from the set of coloured graphs with vertex set V to the set of graphs with 
vertex set V such that, for any coloured graph (G, tt), £(G, tt) = t(G) for some permutation 
r with t(7t) = c(ir) and £(o~(G), <x(tt)) = £(Cr, 7r) for every permutation a of V. 

Theorem. [M, Theorem 2.2], [N, Theorem 1] Let {G,ix) and (H,p) be coloured graphs on 
the common finite vertex set V such that n and p have the same number of vertices in the i th 
colour class, for each i. Further, let £ be a canonical labelling map for graphs with vertex set 
V. Then £(G, tt) = £(H, p) if and only if o~(G) = H for some colour-preserving permutation 
a. 

Initially, nauty colours a graph with a single colour and then refines this colouring using 
a specified vertex invariant. Then, using the above theorem, nauty can determine if two 
graphs are isomorphic by checking if their canonical labellings, after refinement, are the 
same. The power of nauty is that it implements fifteen different vertex invariants (see [N, 
Section 9]), each of which is more or less useful depending on the class of graphs being tested. 

Examples. We first load Nauty, which automatically loads Edgeldeals for access to the 
Graph class. The nauty software stores graphs in two different string formats, Graph6 and 
Sparse6 (see [N, Section 19] for a complete description); both formats are handled by Nauty 
in the method stringToGraph. The method graphToString always returns the Graph6 
string- representation of the graph; e.g., the five-cycle can be represented as "Dhc" and the 
complete graph on five vertices is represented as "D~{" in the Graph6 string-representation. 

11 : needsPackage "Nauty"; 

12 : R = QQ[a. . e] ; 

13 : graphToString cycle R 

03 = Dhc 

14 : graphToString completeGraph R 

04 = D~{ 

15 : edges stringToGraph (" Dhc " , R) 

05 = -[{a, b}, {b, c}, {c, d}, {a, e}, {d, e}} 

Caveat. Due to the way graphs are constructed in Edgeldeals, the conversion process from 
string format to graph format is costly. Thus, when many consecutive manipulations will 
be done through Nauty, we recommended that the graphs be left in string format until it is 
necessary to convert back to graph format. 

The most powerful feature of Nauty is determining if two graphs are isomorphic. We can 
demonstrate it for two given graphs with arelsomorphic; further we can reduce a list of 
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graphs to be pairwise non-isomorphic with removelsomorphs. We demonstrate the latter 
on a list of 5! = 120 different labellings of the five-cycle, which is represented as "Dhc" in 
the Graph6 string-representation. 

16 : G = graph {{a, c}, {c, e}, {e, b}, {b, d}, {d, a}}; 

17 : arelsomorphic (cycle R, G) 

07 = true 

18 : removelsomorphs apply (permutations gens R, P -> 

graphToString graph apply(5, i-> {P_i 3 P_((i+l)°/.5)})) 

08 = -[Dhc} 

Nauty also includes methods for generating all graphs on a given number of vertices, 
possibly with restrictions to simple properties such as the number of edges, with the method 
generateGraphs. We generate all graphs with between one and nine vertices and verify the 
counts with the Online Encyclopedia of Integer Sequences [OEIS, A000088]. 

Caveat. As nauty does not handle graphs with zero vertices, Nauty will throw an error 
when graphs with zero vertices are requested or encountered. 

19 : A000088 = apply(1..9, n -> #generateGraphs n) 

09 = (1, 2, 4, 11, 34, 156, 1044, 12346, 274668) 

We can also select all graphs in a list with a given property. To do this, we first build a 
filter with the method buildGraphFilter and then use the method f ilterGraphs to select 
all graphs in the list which pass the filter; see the documentation for buildGraphFilter to 
see all the possible properties which can be filtered for. First we generate all bipartite graphs 
with between one and twelve vertices. We then filter this list for forests, i.e., graphs without 
cycles. We verify the counts with [OEIS, A005195]. 

110 : B = apply(1..12, n -> generateGraphs (n, OnlyBipartite => true)); 

111 : forestsOnly = buildGraphFilter {"NumCycles" => 0}; 

112 : A005195 = apply(B, graphs -> #f ilterGraphs (graphs , forestsOnly)) 
ol2 = (1, 2, 3, 6, 10, 20, 37, 76, 153, 329, 710, 1601) 

Caveat. When filtering graphs for connectivity, Nauty uses the following definition of k- 
connectedness. Specifically, a 0-connected graph is a disconnected graph and, for k > 0, a 
/c-connected graph is a graph that can be disconnected by removing k vertices but not by 
removing k — 1 vertices. 

Thus, if we wish to select only connected graphs, then we must create a filter for graphs 
which are not 0-connected. Alternatively, the method generateGraphs has the option 
OnlyConnected which forces only connected graphs to be returned. We demonstrate the 
former by filtering the list of bipartite graphs for trees, i.e., connected graphs without cycles. 
We verify the counts with [OEIS, A000055]. 

113 : treesOnly = buildGraphFilter {"NumCycles" => 0, 

"Connectivity" => 0, "NegateConnectivity" => true}; 

114 : A000055 = apply (B, graphs -> #f ilterGraphs (graphs , treesOnly)) 
ol4 = (1, 1, 1, 2, 3, 6, 11, 23, 47, 106, 235, 551) 

Last, when testing conjectures — especially when dealing with graphs on many vertices — 
it is nice to be able to generate large lists of random graphs. Edgeldeals provides the 
method randomGraph, but this method only generates graphs with a specified number of 
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edges and, moreover, generates only one graph at a time. Nauty provides the method 
generateRandomGraphs, which randomly generates multiple graphs on a fixed number of 
vertices with a specified edge probability 

Erdos and Renyi showed in [ER, Theorem 1] that a random graph on n vertices with 
edge probability ( 1+ ^ logra is almost always connected while a graph on n vertices with edge 
probability ^ 1 ~ £ ^ logn is almost never connected, at least as n tends to infinity. We demonstrate 
this property by generating 100 random graphs on n vertices where e is 1 in the first case 
and \ in the second. 

115 : connected = buildGraphFilter -["Connectivity" => 0, 

"NegateConnectivity" => true}; 

116 : prob = n -> log(n)/n; 

117 : apply(2..30, n-> #f ilterGraphs ( 

generateRandomGraphs (n, 100, 2* (prob n)), connected)) 

017 = (70, 83, 88, 96, 94, 97, 97, 96, 97, 96, 99, 96, 99, 98, 98, 

96, 99, 96, 99, 96, 98, 97, 97, 98, 100, 99, 97, 99, 99) 

118 : apply(2..30, n-> #f ilterGraphs ( 

generateRandomGraphs (n, 100, (prob n)/2), connected)) 

018 = (17, 8, 3, 2, 3, 0, 1, 0, 1, 2, 3, 0, 1, 0, 0, 1, 1, 2, 0, 

2, 0, 1, 0, 0, 0, 0, 0, 0, 0) 
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